Large Negative Thermal Expansion in Pentacene due to Steric Hindrance 
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The uniaxial negative thermal expansion in pentacene crystals along a is a particularity in the 
series of the oligoacenes, and exeptionally large for a crystalline solid. Full x-ray structure analysis 
from 120 K to 413 K reveals that the dominant thermal motion is a libration of the rigid molecules 
about their long axes, modifying the intermolecular angle which describes the herringbone packing 
within the layers. This herringbone angle increases with temperature (by 0.3-0.6° per 100 K), and 
causes an anisotropic rearrangement of the molecules within the layers, i.e. an expansion in the 
b direction, and a distinct contraction along a. Additionally, a larger herringbone angle improves 
the cofacial overlap between adjacent, parallel molecules, and thus enhances the attractive van der 
Waals forces. 

PACS numbers: 65.40.De, 72.80.Le 



I. INTRODUCTION 

The thermal expansion observed in most solids is a 
hallmark of anharmonic interatomic or intermolecular 
potentials. Whereas a hard, steep potential results in 
relatively small thermal expansion, a softer, shallow po- 
tential leads to substantial thermal expansion, e.g. in or- 
ganic molecular materials, where the molecules are held 
together by weak van der Waals (vdW) forces. Thus 
the thermal expansion is caused by changes of the ar- 
rangement of the molecules rather than by the small 
temperature-dependence of the intramolecular (covalent) 
bonds. 

On the other hand, negative thermal expansion (i.e. 
contraction) can not be explained by a simple anhar- 
monic two-center potential. Quite a few inorganic ma- 
terials show either volume or uniaxial negative thermal 
expansion (NTE)i with ZrW208 as the most prominent 
example of large, isotropic NTE over a wide temperature 
rango^i^ii. In Q;-ZrW208 (low temperature phase) and /3- 
ZrW208 (high temperature phase), the linear expansion 
coefficient a is -8.7 x 10"^ K"! and -4.9 x lO'^K'^ 
rcspectivelji^. In contrast, only very few examples of 
(uniaxial) NTE in organic systems are known so far (see 
Ref. and refs. therein). 

In general, NTE is caused by special geometrical ar- 
rangements of the atoms/molecules and peculiarities of 
the bonding which restrict the thermal movement of the 
atoms/molecules such that NTE occurs for geometrical 
reasons. Often, such materials are characterized by open, 
underconstrained lattices or by a network of rigid, un- 
derconstrained entities. For example, rotation of the 
rigid WO4 building blocks causes the overall contrac- 
tion in ZrW208i. For organic materials, different mech- 
anisms may be involved. For instance, polyethylenoS, 



contracts along the covalent chains as a consequence 
of increased vibrations perpendicular to the chains, i.e. 
due to motions in direction of weak vdW interactions^ii^. 
Another mechanism has been observed for a dipeptide 
(TrpGly-Il20) crystal^: The ordering of water molecules 
(located within the supramolecular peptide helices) with 
decreasing temperatures is responsible for the uniaxial 
NTE. 

The intermolecular potential and its anharmonicity de- 
termine both thermal expansion and the details of ther- 
mal libration and translation motions. Both of them 
are accessible to experimental probes. To this end, 
we have performed full x-ray structure analyses of pen- 
tacene single crystals in the temperature range from 120- 
400 K. The thermal expansion ellipsoid is calculated for 
pentacene and the other acene crystals, and is highly 
anisotropic. Pentacene crystals contract upon warming 
along one direction which is close to the in-plane a direc- 
tion. The microscopic origin of NTE is elucidated by an 
analysis of the refined atomic displacement parameters 
in terms of molecular motions. Anomalously large libra- 
tions about the long molecular axis at elevated tempera- 
tures are observed, while the other librations and trans- 
lations are mostly suppressed due to steric hindrance 
within the dense herringbone packed layers. As a con- 
sequence, the herringbone angle changes, resulting in an 
expansion approximately along b and a contraction along 



II. EXPERIMENTAL 

Single crystals of pentacene were grown in a tempera- 
ture gradient in vacuo, or by physical vapor phase trans- 
port with high purity argon as the transport gasiSiii. 



2 




300 320 340 360 380 400 420 
T[K] 



FIG. 1: Temperature dependence of the lattice constants of 
pentacene (single crystal, CuKc-data). Expansion in the ab- 
plane is strongly anisotropic and negative along a. (This data 
corresponds to the open symbols in Fig. (2]) 



Pentacene crystals were investigated between room tem- 
perature and 400 K, using an Enraf/Nonius CAD4 
diffractometer (CuKq) with a home-made temperature 
control system, and additionally between 120 K and 
413 K on the Swiss-Norwegian Beamline (SNBL, A = 
0.71 A) at ESRF (Grenoble). 

Full structure data for pentacene single crystals at dif- 
ferent temperatures are available from the CCDC with 
ref.no. xxxxxx-yyyyyy or upon request from the au- 
thors. Crystallographic data for naphthalene, anthracene 
and tetracene were taken from the literatur o^^'^'^i^^'^^ . 

The thermal expansion was calculated from the change 
of the unit cell parameters with temperature. As the var- 
ious acenes crystallize in different, non-orthogonal unit 
cells, a direct comparison of the unit cell parameters is 
of limited value. More appropriate is a transformation of 
the thermal expansion along the crystal axes (a^) into an 
orthogonal reference system, yielding the direction and 
the magnitude of the thermal expansion ellipsoid main 
axes a'^ (principal axes, eigenvalues of the thermal expan- 
sion tensor), following a procedure described earlie r^^i^^ . 
For the thermal expansion analysis, the CAD4 data are 
considered. 

Additionally, molecular libration and translation pa- 
rameters for pentacene were obtained by carrying out 
a TLS analysisiSii^ from the full structure data (syn- 
chrotron data, refined with anisotropic atomic displace- 
ment parameters), using PLATOr*22,. This analysis is car- 
ried out in the molecular inertial system and assumes 
rigid molecules. Biirgi et al. showed the amplitudes 
of the intramolecular modes to be small compared with 
the rigid body motioi*^. Additionally, molecular dynam- 
ics studies indicate that in the energy range of interest 
the excited modes are to a high degree of intermolecular 
type2^. 



FIG. 2: The eigenvalues of the expansion tensor (a^) for 
crystals of acenes with n— 2, 3, 4, and 5 benzene rings. In 
pentacene the expansion coefficient Q3 is negative. The eigen- 
vectors are not parallel to the crystal axes (except for one 
direction along the b axis in the monoclinic systems.) The 
two datasets for pentacene (■, □) are both measured with 
CuKa and point detector, but with different temperature con- 
trols. Data for naphthalene taken from Oddershede et alM- 
and Brock et aL— (light symbols), for anthracene from Brock 
et alt^ and for tetracene from Ref. [Tsl . 



III. RESULTS AND DISCUSSION 

For all investigated crystals, a structure identical with 
the '"H"' polymorp h^^i^^i^^ was observed. A plot of the 
unit cell axes in pentacene reveals a negative thermal ex- 
pansion coefficient along a (Fig. [1]). A contraction 
of the unit cell along a can also be seen in the pow- 
der data of Mattheus et al^, although the single crystal 
data in Ref. [2^ would not indicate NTE. The thermal ex- 
pansion coefficients (expressed in orthogonal principal 
axesi^iii) are plotted in Fig. [2] for two pentacene sam- 
ples and the shorter acenes (Refs. |Il|lilI3lil) . Most 
outstanding is the behavior of pentacene: the ther- 
mal expansion in one direction is distinctly negative 
(—35 ± 6xlO~^K~^), and positive in the two other di- 
rections. The trend for naphthalene, anthracene and 
tetracene indicates a dependence of the largest thermal 
expansion coefficient on the length of the molecule. There 
is at least one direction with a value close to zero in all 
three. Pentacene is anomalous among the series not only 
in terms of the NTE, but also when its relatively large 
volume expansion coefficient (=X]'^D considered. 

Fig. [3] shows the direction of the principal axes of the 
thermal expansion ellipsoid with respect to the pentacene 
molecules and the unit cell. The direction of most pro- 
nounced NTE (axis 3, 013) encloses an angle of 22° with 
the crystallographic axis a and is nearly perpendicular 
to the crystallographic axes b and c. The direction of 
maximal expansion (axis 1, a[) encloses an angle of 26° 
with b and 50° with c, being almost perpendicular to a. 
The anomalously large intermediate eigenvalue (axis 2, 
a'2) corresponds to a direction out of the a6-plane, but is 
not exactly parallel to the long molecule axis. 
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FIG. 3: Thermal expansion in pentacene: The main axes of 
the thermal expansion ellipsoid are non-parallel to all crystal 
axes. Axis 3 (eigenvector to the value of minimal expansion) 
corresponds to a negative expansion, i.e. contraction in that 
direction nearly parallel to a, whereas maximum expansion 
directs about along b and towards the center molecule (axis 
1). For numerical values, see Fig. [2] ^ is the herringbone 
angle, and 5 the tilt angle. 



To explore the microscopic cause of the negative ther- 
mal expansion in pentacene, we further analyze the tem- 
perature dependence of the thermal motions and the ori- 
entation of the molecules (herringbone/tilt angle, defini- 
tions see Fig. [3]). 

The libration and translation parameters for the rigid 
molecule motion were calculated from synchrotron data 
and are listed in Table |T1 The longer acene molecules 
show the expected increase in the anisotropy of the 
molecules' inertial tensor. The longer the molecule the 
more suppressed ( "frustrated" ) are the librations due to 
steric hindrance, and the more closely docs Ln coincide 
with the molecular inertial axis 1 (given by /i ) . In light of 
the negligibly small off-diagonal elements of the libration 
tensor (given by the coupling between the motions), Ln 
describes a movement independent from ^22/331 in 
is aligned with Ii only in the case of pentacene (misalign- 
ment < 5°) and tetracene (<7°). The translations T^, 
however, remain almost constant with increasing num- 
ber (n) of benzene rings (Tab. |T|. 

The herringbone angle 9 (see Fig. [3] for definition) is a 
central parameter to describe the packing within the lay- 
ers, and is essentially the same for all acenes (51-53° at 
room tempcratur), irrespective of the symmetry of the 
unit cell, and shows the same qualitative temperature 
dependence. With increasing temperature, the herring- 
bone angle increases by 0.3-1° per 100 K, with a more 
pronounced change in naphthalene and anthracene than 
in tetracene and pentacene. 

The increase of the herringbone angle and the increas- 
ingly dominant librations Ln about the long molecular 
axis suggest a distinctly anisotropic expansion in the ab- 
plane, as sketched in Fig. [51 The librating molecules 
at the corners of the unit cell interact with the one in 
the center, pushing each other further apart in h direc- 



TABLE I: Translation and libration parameters for linear 
acenes at room temperature, calculated in the molecules' in- 
ertial axes coordinate system. Libration parameters La are 
in deg^, translation parameters Ti in 10~^ A^. For the mono- 
clinic acenes, the two molecules in the unit cell are crystallo- 
graphically identical. In case of triclinic symmetry, they are 
independent, creating additional parameters L^j and T/ for 
the second molecule. The main axes of the inertial tensor, 
are given in A^-amu, and also normalized with the number n 
of benzene rings. Ji (smallest momentum) points along the 
long molecule axis. 



(n) 


Naph" (2) 


Ac'' (3) 


Tc" (4) 


Pc (5) 


Symmetry 


mono 


mono 


tri 


tri 


h (-n) 
h i^n) 
h i^n) 


130 (65) 
356 (178) 
486 (243) 


189 (63) 
1001 (334) 
1190 (397) 


252 (63) 
2154 (539) 
2406 (602) 


311 (62) 
3957 (791) 
4268 (854) 


Lu 


19.32 


15.3 


11.07 


9.55 


L22 


12.18 


7.2 


4.51 


2.69 




16.55 


10.3 


3.71 


2.17 


i'll 






13.01 


9.59 


T' 

^22 






5.25 


2.76 


Li3 






4.19 


2.06 


Tl 


3.18 


4.77 


5.461 


4.240 


T2 


2.37 


3.02 


2.991 


2.288 


T3 


1.58 


1.94 


2.632 


2.205 


Ti 






4.636 


4.504 


n 






2.876 


2.408 


n 






2.646 


2.069 



"Values from Ref,i^. 
''Values from Ref™. 
"^Values from Ref™. 
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FIG. 4: Temperature dependence of the molecular librations 
in pentacene. The libration about the long molecule axis, 
ill, is most pronounced, and strongly increases with temper- 
ature. The mean rotational amplitude is approximately ^/La, 
i.e. > 3° about the long axis at room temperature. Inset: The 
measured (static) herringbone angle 9 (•) and the mean ther- 
mal libration angle ±\/iii (connected by the broken line) as 
a function of temperature. If the thermal expansion along b 
was caused by a static rotation only, the corresponding angle 
(o) varied more rapidly than the measured herringbone angle 
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tion as the libration amplitude grows with temperature. 
Thus, steric hindrance between the herringbone-packed 
molecules is central to the NTE in pentacene. 

It is noteworthy that this simplified picture of one dom- 
inant libration in is most appropriate for pentacene and 
tetracene, where the inertial tensor is anisotropic enough 
to align ill with the long molecular axis within 5-7°. 
Lattice dynamics calculations indicate that the lowest en- 
ergy modes involve essentially pure rotational and trans- 
lational iniermolecular motions, whereas miramolecular 
excitations require higher energy (e.g. Tab. II in Ref. 
[22I ). Thus Fig. [5] captures the molecular motions most 
relevant for thermal expansion in pentacene. 

The observed difference of a more anisotropic (and 
even negative) thermal expansion in pentacene compared 
to tetracene can be qualitatively understood as the result 
of two counteracting types of thermal motion: in-plane 
translations and the librations along the long molecular 
axis. The libration apparently dominates in pentacene. 

Assuming an effective width of the molecules of 7 A, 
the thermal expansion along 6, A6caic, originating from 
the change of the static (i.e. average) herringbone an- 
gle can be easily calculated geometrically, as suggested 
in Fig. [5l (The tilt of the molecules with respect to b has 
been taken into account.) Since A6caic < A^cxp, there 
is obviously an additional contribution from dynamic ef- 
fects, i.e. the librations are partially out-of-phase. In a 
different picture, this corresponds to lattice phonons with 
q 7^ 0. Consequently, an effective herringbone angle ^off 
was calculated, which varies more rapidly with temper- 
ature than the measured, "averaged" herringbone angle 
9, as shown in the inset of Fig. |4l Taken together, the 
information in Fig. [4] reflects the anharmonic potential 
sensed by the thermal motion. 

An additional factor for the negative expansion is the 
increased vdW force along a due to an increased cofa- 
cial overlap of the pentacene molecules: The tilt angle 5 
against c* (definition in Fig.[3|) decreases from 6=25.4° at 
120 K over 6=24.7° at 293 K to 5=24.0° at 413 K (aver- 
aged for the two inequivalent, quasi-parallel molecules). 
As a consequence, the displacement of the pentacene 
molecules is slightly reduced from 2.10 A at 120 K over 
2.05 A at 295 K to 1.98 A at 413 K, and thus the cofacial 
overlap and the attractive vdW forces along a increase. 

It can be expected that the asymmetry and the uniax- 
ial negative expansion is even more pronounced in longer 
herringbone packed acenes, e.g. hexacenc. In particu- 
lar the effect of increasing vdW forces along the short 
unit cell axis, and the related change of the tilt angle 



are enhanced, as well as the steric hindrance within the 
herringbone layer. However, the decreasing stability and 
the increasing ease of oxidation with increasing n, as well 
as the reported difficulties in growing hexacene crystals 
with reasonable quality^ make single crystal studies a 
challenge. 




FIG. 5: Scfiematic view of the anisotropic expansion in the 
a,b plane. The libration along the long molecular axis, Ln, 
together with a change of the herringbone angle, are responsi- 
ble for the anisotropic thermal expansion, with a contracting 
and b expanding. The other thermal libration and transla- 
tion motions are strongly suppressed ("frustrated") by steric 
hindrance. 



IV. CONCLUSIONS 

In summary, the uniaxial negative thermal expansion 
in pentacene along a is a particularity in the series of the 
acenes. Pentacene contracts along a as a result of the 
increasing herringbone angle (A0=O.3-O.6° per 100 K), 
supported by the frustation of larger librations and trans- 
lations perpendicular to the long molecule axis due to 
steric hindrance. Furthermore, a reduction of the relative 
shift of the molecules along a with temperature increases 
the attractive vdW forces along a. 
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